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We present an algorithm that provides a Grobner representation of a binary linear code 
C and also the set of all coset leaders, denoted by CL(C). Its complexity is linear on the 
number of elements of CL(C), which is smaller than exhaustive search in FS. It is also 
shown that the same algorithm could provide a set of codewords which form a test set 
for decoding by a gradient-like decoding algorithm. 
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1. Introduction 

Complete decoding problem (CDP) addresses the question of retrieving a codeword 
that minimized the hamming distance with respect to a given received word. On the 
other hand the t-bounded distance decoding problem (t-BDD) looks for a codeword, 
if possible, whose hamming distance to a given received word is at most t. Observe 
that if t coincides with the covering radius of the code C, then both problems over C 
are equivalent. In fact in the CDP of a linear code, those errors that can be corrected 
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are those related to the coset leaders since if there is more than one leader in a coset 
there is more than one choice for the error. 

Let H be a parity check matrix of a linear code and s be the syndrome of a 
received word. It is widely known that t-BDD problem is equivalent to the coset 
weights problem (CWP) which can be stated as the problem of finding an error 
vector of hamming weight at most t with the same syndrome as the received word. 
If we express the problem schematically it becomes: 

Input: A binary r x n matrix H, a vector s G F£ and a non-negative integer t. 
Problem: Does a binary vector e £ F£ of Hamming weight at most t exist such 
that He = s? 

Therefore finding a solution the previous problem is equivalent to solving the 
£-BDD problem of the received word. 

The set of coset leaders in linear codes has been related to the set of minimal 
support codewords which have been used in maximum likelihood decoding analysis 
|2I14) and in secret sharing schemes since they describe the minimal access structure 
|15j . Moreover a slight modification of the main algorithm proposed in this paper 
to compute the set of all coset leaders can be done to obtain a test set (called leader 
codewords) that allows gradient descent decoding. In [5] the authors proved that the 
two gradient descent decoding algorithm for binary codes proposed independently 
by Licbler [13) and by Ashikhmin and Barg [1] are equivalent. 

All problems mentioned above are considered to be hard computational prob- 
lems (see for example [213] ) even if preprocessing is allowed [9]. Our purpose is to 
extend the previous work on Grobner representation to obtain the set of all coset 
leaders taking advantage of the additive structure of the cosets which could be done 
efficiently. 

Note that the structure described in this paper is related to the monotone struc- 
ture of the sets of correctable and uncorrectable errors and also introduced in [11] . 
They introduce trial sets of codewords, considering the linear ordering of all vectors 
of a fixed length defined by x ^ y if and only if x has smaller Hamming weight than 
y or they have the same Hamming weight but is lexicographically smaller than y. 
In that paper they describe the minimal uncorrectable errors under the ordering 
and the so-called larger halves of minimal codewords as well as use the description 
to give a gradient-like decoding algorithm. The same approach is considered in [16] . 

The following is an outline of the paper. In Section [2] we have compiled some 
basic facts on coding theory and the Grobner representation of binary linear codes. 
As for prerequisites, the reader is expected to be familiar with these topics. However 
we will touch only an a few aspect of the theory of Grobner bases since the paper 
is written in a "Grobner bases" -free context. For a deeper discussion of Grobner 
representation for linear codes we refer the reader to |6I8I4| where recent results and 
some applications are indicated, in order to get a general picture on the subjects 
we recommend [TJ. 

In Section [3] our main result is stated and proved. Moreover the main theorem 
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provides an algorithm to compute the set of all coset leaders. The example presented 
at the end of this section of a binary linear code with 64 cosets and 118 coset 
leaders suggests extra applications of the algorithm such as how to obtain the weight 
distribution of the coset leaders or the Newton and Covering radius of a code. These 
applications do not pose a large additional cost to the proposed algorithm. 

Section |4] is devoted to show how the previous algorithm can be adapted to 
compute a test set for the code which we refer to as leader codewords. Not only 
do we prove that they are zeroneighbors but also that the knowledge of the set of 
leader codewords can be used to compute all coset leaders corresponding to a given 
received word. 

In Section [5] some complexity results are discussed. Finally some conclusions and 
new directions for future research are outlined. 

2. Preliminaries 

By F2 we denote the finite field with 2 elements. A linear code C over F2 of length 
n and dimension k, or an [n, k] binary code for short, is a fc-dimcnsional subspacc 
of FJ. We will call the vectors v in FJ words and the particular case where v G C, 
codewords. For every vector y £ its support is define as its support as a vector 
in FJ, i.e. supp(y) = {i \ y t ^ 0} and its Ramming weight, denoted by wjj(y) as 
the cardinality of supp(y). 

The Hamming distance, dn (x, y), between two vectors x, yeFJ is the number 
of places where they differ, or equivalently, <if/(x, y) = w#(x — y). The minimum 
distance d of a linear code C is defined as the minimum weight among all nonzero 
codewords. 

Definition 1. The words of minimal Hamming weight in the cosets o/FJ/C are 
the set of coset leaders for C inFJ. We will denote by CL(C) the set of coset leaders 
of the code C and by CL(y) the subset of coset leaders corresponding to the coset 
C + y. 

Remark 2. It is well known that CDP over the code C has a unique solution for 
those vectors in the union of the Hamming balls of radius t around the codewords 
of C i.e. B(C,t) = {y G FJ | c G C s.t. d H (c,y) < t}, where t = [^J is the error- 
correcting capacity of C and ['J denotes the greatest integer function. 

From now on {e; | i = 1, . . . ,n} represents the canonical basis of FJ . The fol- 
lowing theorem gives us a nice relationship between the coset leaders. 

Theorem 3. Let w G CL(C) such that w = y + e^ for some word y G FJ and 
i G supp(w), then y G CL(w). 

Proof. See [H Corollary 11.7.7]. □ 
Definition 4. A Grobner representation ofW^/C is a pair (N,(/)) where: 
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(1) N is a transversal for the cosets of FJ/C such that G N and for each n G 
-/V \ {0} there exists i G {1, . . . ,n} such that n = n' + e^ with n' G N and 

(2) (f> : N x {e;}" =1 — > N is a function called Matphi that maps each pair (n,e.;) 
to the element of N representing the coset o/n + e^. 

3. Computing the set of coset leaders 

We will use the following characteristic crossing functions: 

T : N s — ► F§ and A : F| — ► N s . 

The maps T is reduction modulo 2 whereas the map ▲ replaces the class of 0, 1 
by the same symbols regarded as natural numbers. The integer s is determined by 
context and both maps act coordinate-wise. 

Definition 5. An ordering -< on F£ is a weight compatible ordering if for any 
a, b £ we say a -< b if 

w#(a) < w#(b) or w#(a) = w#(b) and Aa -<i Ab 

where -<i is any admissible order on N n . 

Remark 6. A weight compatible ordering -< is in general not admissible on Fj but 
it satisfies the following conditions: 

(1) -< is a Noetherian order on Fj since Fj is finite. 

(2) Let v, w G ¥2 , if supp(v) C supp(w) then v -< w. 

Definition 7. The object List is an ordered set of elements in Fjj 1 w.r.t. a weight 
compatible order -< verifying the following properties: 

(1) G List. 

(^j If v E List and wh-(v) = wjj (iV(v)) t/ien {v + e^ i ^ supp(v)} C List, 
where N(w) = min^ {w | w G List PI (C + v)}. 

Remark 8. If the second condition holds for v G F£ then v is a coset leader for 
the coset C + v. 

Definition 9. For eac/i v G List we have defined N(v) as the minimal element 
w.r.t. a fixed weight compatible order -< which belongs to the intersection List PI 
(C + v). FFi/Z denote by AT the set of distinct N(v) with v G List. 

Next theorem states that the object List includes the set of coset leaders of a 
given binary linear code. 

Theorem 10. Let w G F§. I/w G CL(C) f/ien w G List. 

Proof. We will proceed by Noetherian induction on FI? with a weight compatible 
ordering -<. 
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The statement is true for S Fj by definition. Now assume that the desired 
property is true for any word u smaller than an arbitrary but fixed w G CL(C) \ {0} 
w.r.t. -<, i.e. 

if u G CL(C) and u -< w then u G List. 

Let i G supp(w), we can write w = u + e; with i ^ supp(u), or equivalently, 
supp(u) C supp(w), hence u -< w. In addition, since w G CL(C), then by Theorem 
[3] also u G CL(C). So if we invoke the induction hypothesis we have that u G List. 
Moreover w#(u) = wg (iV(u)) which is clear from the fact that u G CL(C). We now 
apply property 2 of Definition [7] which gives, as claimed, that w = u + G List 
with i £ supp(u). □ 



Algorithm 1: CLBC Algorithm 



0; CL(C) <— 0; S 



tH J 



t = t' + e fe with t'ejV 



Data: A weight compatible ordering -< and a parity check matrix H of a 
binary code C. 

Result: The set of coset leaders CL(C) and (J\f, <j>) a Grobncr representation 
for C. 

1 List i — [0]; Af 

2 while List ^ do 

3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 
16 



t < — NextTerm[ListJ 
j < — Member[s, <S]; 
if j 7^ false then 
for k G supp(t) : 
| ct>(t',e k ) <— t 
endfor 

if wjj(t) = wn{tj) then 

CL(C)[tj] ^— CLC[tj] U {t}; 
List ^ — InsertNext[t, List]; 
end if 
else 

r i — r + 1; t r < — t; Af < — jVU {t r }; 
CL(C)[t r ] < — {t r }; 5 < — S U {s}; 
List = InsertNext[t r , List]; 
for fc G supp(t r ) : t r = t' + e^ with t' G Af 
0(t',e fc ) < — t r ; 
</)(t r ,e fe ) < — t'; 
endfor 
end if 



17 
18 
19 
20 
21 

22 end while 



The subfunctions used in the algorithm are: 
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• InsertNext[t, List], adds to List all the sums t + with k ^ supp(t) and 
keeps List in increasing order w.r.t. the ordering -<. 

• NextTerm[List], returns the first element from List and deletes it. If List is 
empty returns 0. 

• Member[obj, G], returns the position j of obj in G if obj G G and false other- 
wise. 

Remark 11. First we perform subroutine t = NextTerm[List] where the element 
t is deleted from the set List. Then subroutine InsertNext[t, List] is carried out 
which inserts in List all the elements of the form: 

t' = t + efe with k £ supp(t), i.e. t' y t. 

Therefore all the new elements inserted in List are greater than those that have 
already been deleted from it with respect to -<. 

Theorem 12. The CLBC algorithm computes the set of coset leaders of a given 
binary code C and its corresponding Matphi function. 

Proof. Let us first prove that the set List is well defined according to Definition 
[7l G List verifying property 1, by Step 1. In Step 3 the syndrome of t = 
NextTerm[List] is computed, then we have two possible cases based on the outcome 
of Step 4: 

(1) If j = false then the coset C + t has not yet been considered. Thus, according 
to Remark [TT1 we have that N(t) = t and Step 16 guarantees property 2. 

(2) Contrarily, if j ^ false, then the element N(t) = tj has already been com- 
puted. However, if t G CL(C), or equivalently, w#(t) = w#(tj), then Step 11 
certified property 2. 

Therefore Algorithm CLBC construct the object List in accordance with Definition 

m 

Furthermore, on one hand Step 10 and 15 assure the computation of the set of 
coset leaders of the given code; and on the other Step 7, 18 and 19 computes the 
Matphi function. Note that Step 17 is necessary since the above first case assure 
that N(t r ) = t r so by Theorem [3] t' G CL(C). But t r = t' + e k with k G supp(t r ) 
so by Remark [IT] t' -< t r has already been considered on the algorithm. Thus we 
have actually proved that the CLBC algorithm guarantees the desired outputs. 

Finally notice that the cardinality of the set List is bounded by n times the 
cardinality of CL(C) and Step 11 and 16 guarantee that when the complete set 
of coset leaders is computed no more elements are inserted in List while Step 3 
continues deleting elements from it. Thus after a finite number of steps the set List 
get empty. Consequently, Step 2 give the end of the algorithm. □ 

Remark 13. Note that the Algorithm CLBC returns (A/", (f>) that fulfill Definition 
@J for correctness we refer the reader to [5J Theorem 1] . Furthermore, by Definition 
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[9] those representative of the cosets given by M are the smallest terms in List w.r.t. 
-<. 



Remark 14. The collection of programs and procedures GBLA_LC (Grobner 
Basis by Linear Algebra and Linear Codes) has been already presented in previous 
works (sec for instance |6I8I4| ). This framework consist of various files written in 
the GAP [TI5] language and included in GAP's package GUAVA 3.10. 

Now we have implemented the CLBC algorithm and added to the collection 
GBLAJLC. This function GBLA. has the same input of the algorithm and returns 
a list with three components given by the set of coset leaders with respect to any 
total degree compatible ordering, the function Matphi and the error correction 
capacity of a given binary code. 

Example 15. Consider the [n = 10, k = 4, d = 4] binary code C defined by the 
following parity check matrix: 

/ioooiooooo\ 

10 110 10 



110 10 10 
1110 10 
11110 10 

yi liioooooiy 



T 6xl0 



The CLBC function returns the following list of cosets leaders with 64 compo- 
nents where each components corresponds to a coset with its coset leaders ordered 
w.r.t. the degree reverse lexicographic order -<. Thus the set A/" is composed by the 
first element of each component. Recall that this function also returns the function 
Matphi and the error correction capacity of C, in this case t = 1. We denote by 
CL(C)^ the j-th element of the set of coset leaders of weight i. 

The main difference between this paper and previous works is the consideration 
of all coset leaders and not just those belonging to J\f. Note that no subword of 



two elements of y = e4 + e§ + eg G 



CL(C)§ 



is part of A/", i.e. e± + e 5 G CL(C)2, 



e4 + eg G CL(C)5> 3 and e$ + eg G CL(C);j, do not lie in Af. Therefore the importance 
of the second property of the Definition [7] to obtain the complete set of coset leaders. 

The algorithm could be adapted without incrementing the complexity to get 
more information such as: 

• The Newton radius v(C) of a binary code C is the largest weight of any vector 
that can be uniquely corrected, or equivalcntly, v(C) is the largest value among 
the cosets with only one coset leader. In our example it suffice to analyze the last 
element of the list CL(C) to obtain the coset of highest weight which contains 
only one leader, i.e. v(C) = 3 since CL(C)g 3 = [e 5 + e 9 + ei ]. 

• The Covering radius p{C) of a binary code C is the smallest integer s such that 
¥2 is the union of the spheres of radius s centered at the codewords of C, i.e. 
p(C) = maxygpn min cS c rfff(y, c). It is well known that p(C) is the weight of 
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h 



Coset Leaders CL(C) 


CL(C) 


[0] 


CL(C)i 


feil, 

L -LJ ' 


[e 2 ], [e 3 ], [e 4 ], [e 5 ], [e 6 ], [e 7 ], [e 8 ], [e 9 ], [e 10 ], 






h e 2 ,e 5 + e 6 ], 


[ex 4 e 3 ,e 5 4 e 7 ], [ex 4 e 4 ,e 5 4 e 8 ], 






1- e 5 ,e 2 + e 6 ,e 3 + e 7 , e 4 + e 8 ], [ex + e 6 , e 2 + e 5 ], 






he 7 ,e 3 + e 5 ], 


[ex 4 e 8 ,e 4 4 e 5 ], [ex4e 9 ], [ex4exo], 






h e 3 ,e 6 + e 7 ], 


[e 2 4 e 4 ,e 6 4 e 8 ], [e 2 4 e 7 , e 3 4 e 6 ], 


CL(C) 2 




1- e 8 ,e 4 + e 6 ], 


[e 2 4 e 9 ], [e 2 4 exo], 






he 4 ,e 7 4e 8 ], [e 3 + e 8 , e 4 + e 7 ], [e 3 + e g ], 




fe-? - 


l-exo], 


[e 4 + e g ] , [e 4 + e 10 ] , [e 5 + e g ] , 




fe^ - 


l-e 10 ], [e 6 + e 9 ], [e 6 + e 10 ], [e 7 + e 9 ], 






h exo], 


[es + eg], [e 8 + exo], [eg + exo], 




fei - 


he 2 + e 3 ,ex + e 6 + e 7 ,e 2 + e 5 +e 7 ,e 3 + e 5 + e 6 ], 




[ei - 


Fe 2 + 


e 4 ,e 4 4 


e 6 4e 8 ,e 2 4 e 5 4e 8 ,e 4 4 e 5 4e 6 ], 




[ei- 


he 2 + 


e 7 ,ex 4 


e 3 4e 6 ,e 2 4 e 3 4e 5 ,e 5 4 e 6 4e 7 ], 




[ei- 


Fe 2 + 


e 8 ,ex 4 


e 4 4 e 6 , e 2 4 e 4 4 e 5 ,e 5 4 e 6 4 eg], 




[ei- 


he 2 + 


eg,e 5 4 


e6 4 e 9 ], [ex 4 e 2 4 exo, e 5 4 e 6 4 exo], 




[ei- 


he 3 + 


e 4 ,ex 4 


e 7 4 e 8 , e 3 4 e 5 4 e 8 ,e 4 4 e 5 4 e 7 ], 




[ei- 


Fe 3 + 


e 8 ,ex 4 


e 4 4 e 7 , e 3 4 e 4 4 e 5 , e 5 4 e 7 4 e 8 ], 




[ei- 


he 3 + 


eg,e 5 4 


67 + e 9 ], [ex4e 3 + exo,e 5 + e7 + exo], 


CL(C) 3 


[e x - 


he 4 + 


eg,e 5 4 


es + e 9 ], [ex 4 e 4 4 e 10 , e 5 4 e 8 4 exo], 




[ei- 


he 5 + 


e 9 ,e 2 4 


e 6 4e g ,e 3 4 e 7 4e 9 ,e 4 4 e 8 4 eg], 




[ei - 


l-e 5 + 


exo, e 2 - 


f e 6 4 e 10 , e 3 4 e 7 4 e 10 , e 4 4 e 8 4 exo], 




[ei- 


Fe 6 + 


eg, e 2 4 


eg 4 eg], [ex 4 e 6 + exo,e 2 + e 5 + exo], 




[ei- 


Fe 7 + 


eg,e 3 4 


e 5 4 eg] , [e x 4 e 7 + exo, e 3 + e 5 + exo] , 




[ex- 


Fe 8 + 


e 9 ,e 4 4 


©5+ eg], [ex 4 e 8 + exo,e 4 + e 5 + exo], 




[ex- 


Fe 9 + 


e xo], 






[e 2 - 


he 3 + 


e 8 ,e 2 4 


e 4 4 e 7 , e 3 4 e 4 4 e 6 , e 6 4 e 7 4 eg], 




[e 5 - 


l-e 9 + 


exo] 





the coset of largest weight. Likewise, in our example p{C) = 3 since CL(C)| 3 = 
[eg 4 e 9 4 exo] is the coset of highest weight. 
• The Weight Distribution of the Coset Leaders of a binary code C is the list 
WDCL = (ao, . . . ,a„) where with 1 < i < n is the number of cosets with 
coset leaders of weight i. Note that the set TV is enough to compute this pa- 
rameter. It is clear that 



WDCL = [ 1, 10, 30, 23, 0, 0, 0, 0, 0, 0] 
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• The number of coset leaders in each coset : 



tt(CL) 



"1, 

1,1,1,1,1,1,1,1,1,1, 

2,2,2,4,2,2,2,1,1,2,2,2,2,1,1,2,2,1,1,1,1,1,1,1,1,1,1,1,1,1, 
.4,4,4,4,2,2,4,4,2,2,2,2,4,4,2,2,2,2,2,2,1,4,1 



Note that there are 30 of the 64 cosets where the Complete Decoding Problem 
(CDP) has a unique solution. It is also interesting to note that among the cosets 
with one leaders there are more cosets exceeding the error correction capacity 
(19) than achieving such capacity (11). 

4. Computing a test set 

Definition 16. A test-set T for a given code C is a set of codewords such that every 
word y either lies in the Voronoi region of the all-zero vector, denoted by D(0), or 
there exist t e T such that w#(y — t) < w#(y). 

In this section we show how Algorithm CLBC can be adapted to compute a 
test-set. 

Definition 17. The set of leader codewords of a given code C with parity check 
matrix H is defined as: 



For efficiency reasons we are just interested in a particular case of the above 
object, when supp(ni + ej) R supp(n 2 ) = 0. 

Remark 18. The difference between the algorithms CLBC and CLBC2 are the 
steps from 10 to 12. 

Theorem 19. The CLBC2 algorithm computes the set of coset leaders and the set 
of leader codewords of a given binary code C. 

Proof. By Remark [18] and Theorem [12] it remains to prove that the CLBC2 algo- 
rithm computes the set of leader codewords. We first observe that all the words of 
the set L(C) are leader codewords since t = t' + e$ and tj are in the same coset, in 
particular t,tj E CL(C)[tfc], and by definition we have 

t, tj e CL(C), supp(t) nsupp(tj) = and i supp(t'). 

Furthermore the algorithm ratifies that we are considering all the leader codewords, 
since the elements of List are ordered by -< then, in each loop we study all leader 
codewords of the form ni + e^ + n 2 with ni, n 2 < t with t G CL(C). □ 
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Algorithm 2: CLBC2 Algorithm 



Data: A weight compatible ordering -< and a parity check matrix H of a 
binary code C. 

Result: The set of coset leaders CL(C) and the set of leader codewords L(C) 
for C. 

1 List < — [0]; Af < — 0; r < — 0; CL(C) <— 0; S ^— 0; L(C) <— 0; 

2 while List 7^ do 
3 
4 
5 

6 

7 

8 

9 
10 
11 

12 
13 
14 
15 
16 
17 



t < — NextTerm[List]; s i — tH T ; 
k < — Member [s, <S]; 
if k 7^ false then 

if Wff(t) = wh (tfc) then 

CL(C)[t fe ]^CLC[t fc ]u{t}; 
List ■< — InsertNext[t, List]; 
end if 

for i £ supp(t) : t = t' + e 4 with t' e CL(C) and i ^ supp(t') 
L(C)<— 

L(C) U {t + tj I tj e CL(C)[t fc ] and supp(t) n supp(tj) = 0} 
endfor 
else 

r < — r + 1; t r « — t; Af i — Af U {t r }; 
CL(C)[t r ] <— {t r };5^5U{s}; 
List = InsertNext[t r , List]; 
end if 



is end while 



Definition 20. We define the subset L 1 (C) o/L(C) as 

L 1 (C) = (n 1+ n,+e i , * ^ supp(ni), ni G CL(C), n 2 6 Af 

{ w H (n x + ej) > w ff (n 2 ) and if (n x + e,) = Hn 2 

Following Definition^ Af denotes the set of minimal elements w.r.t. a weight com- 
patible ordering -< of each coset. 

Remark 21. Note that the condition w#(ni + e,) > wh (n 2 ) is imposed just to 
improve the efficiency of computing this set. Therefore L(C) can be rewritten as 

% $. supp(ni), w H (ni + ei) > w H (n 2 ) 
ni, n 2 £ CL(C) and if (ni + ej = if n 2 

Thus the only difference between the sets L 1 (C) and L(C) is that n 2 G Af instead 
of n 2 G CL(C). In other words, to the element n 2 in L 1 (C), not only do we ask for 
belonging to the set of the coset leaders but also to be the least element in its coset 
according to a fix weight compatible ordering -<. 



L(C) = { m + n 2 



Theorem 22. The subset L X (C) o/L(C) is a test-set for C. 
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Proof. Let us consider a word y G F 2 with supp(y) = {i±, . . . ,i m } C {1, . . . ,n) 
such that y ^ CL(C). There must exist an integer 1 < I < w# (y) such that 

Hi := ei + . . . + e*, G CL(C) and m + e ii+1 £ CL(C). 

It is clear that Z is strictly smaller than wjy(y) since y ^ CL(C). 

Now we define n2 = iV(ni +e^ +1 ). Since ni +ei i+1 ^ CL(C). we have wh(ii 2 ) < 
w ff (ni + e i[+1 ). Thus t = ni + n 2 + e, !+1 G L^C). 

Furthermore: 

|supp(t) nsupp(y)| > w ff (m +e i;+1 ) > w H (n 2 ) > |supp(t) n supp(y) 

where y denotes the relative complement of y in F 2 l and, in consequence w#(y — t) < 
w n(y) which completes the proof. □ 

Remark 23. Since L*(C) C L(C) and by Theorem [22] the subset L : (C) form a test 
set for C, then so is the set L(C). 

The following theorem gives a bound for the weight of a leader codeword of a 
given binary code C. 

Theorem 24. Let c G L(C) then w#(c) < 2p(C) + 1 where p{C) is the covering 
radius of C. 

Proof. Let c G CL(C) then there exists ni, n 2 G CL(C) and i ^ supp(ni) such that 
w//(ni + e.;) > w#(n 2 ) and c = ni + e, + n 2 . Applying the definition of covering 
radius we have that w#(rii), w#(n 2 ) < p, thus w# (c) < 2p+ 1. □ 



Algorithm 3: 



Data: A received vector y G F 2 and the set of leader codewords L(C) of a 
binary code C. 

Result: The subset CL(y) of coset leaders corresponding to the coset C + y. 

1 Computes A^(y) by gradientdikc decoding using the test-set L(C) ; 

2 y <— A(y); S <— {y}; L <— L(C); 

3 while there exists ceL : w# (y — c) = w# (y) do 

4 y< — y-c;S< — 5U{y}; 

5 i ■< — L — {c}; 



6 end while 



Theorem 25. Algorithm^ computes the subset CL(y) of coset leaders correspond- 
ing to the coset C + y for a given received vector y G F 2 . 
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Proof. Let us first prove that every z G CL(y) can be rewritten as z = iV(y) — c 
with c G L(C). Let i G supp(z) then z = ni with i ^ supp(ni) then, by Theorem 
El ni G CL(C). Furthermore we have that 

HN(y) = Hz and w^(ni) < w#(z) = wu(N(y)). 

Thus from the definition of leader codewords, N(y) + (ni + e,) G L(C), or equiva- 
lently z = n x + e, = c - N(y) with y G L(C). 

The proof is completed by noting that Theorem [22] guarantees Step 1. □ 



4.1. Leader codewords and zero neighbors 

In this section we will give a brief review of basic concepts from [21 Section 3] and 
thus establish the relation between Zero- Neighbors and leader codewords of a binary 
code C. 

Definition 26. The Voronoi region of a codeword c G C, denoted by D(c), is defined 
as: 

D(c) = {y 6 K I d H (y, c) < d H (y, c') for all c' G C \ {c}} . 

Note that the set of Voronoi regions of a binary code C partition de space F£ , 
however some points of may be contained in several regions. Furthermore the 
Voronoi region of the all-zero codeword D(0) coincides with the set of coset leaders 
of C, i.e. D(0) = CL(C). 

Definition 27. For any subset A C we define X(A) as the set of words at 
hamming distance 1 from A, i.e. 

X(A) = {y G F£ | min {d H (y, a) : a G A} = 1} . 

We define the boundary of A as 5(A) = X(A) U X(W% \ A). 

Definition 28. A nonzero codeword c G C is called Zero-Neighbor if its Voronoi 
region shares a common boundary with the set of coset leaders, i. e. 

5(D(z)) n *(D(0)) ? 0. 

We will denote by Z(C) the set of all Zero-Neighbors of C that is to say: 

Z(C) = {z G C \ {0} : 5(D(z)) D *(D(0)) ^ 0} 

Note that, if z G C \ {0} satisfies that X(D(0)) n D(z) ^ 0, then z G Z(C). 
Furthermore Z(C) is a test-set for C (see for instance [21 Theorem 3.16]). However 
the only property of the set Z(C) that is essential for successful decoding is 

#(D(0))C (J D ( z )- 

ze-Z(C) 
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Thus if we restrict the set Z(C) to a smallest subset verifying the previous property 
we still have a test-set for C. We will denote such subset of Z(C) by Z m i n (C). Note 
that the set iJ m i n (C) may be not unique, however its size is well defined. 

Theorem 29. Let C be a binary code and z G C \ {0}. Then the following are 
equivalent: 

(1) #(D(O))nD(z)^0. 

(2) z 6 L(C). 

Proof. If Af(D(0)) n D(z) ^ then there exists ni G D(0) = CL(C) and i £ 
supp(ni) such that ni + e^ G »Y(D(0)) and ni + e, G D(z). That is to say, 

w#(z - (ni + Bi)) < w H (c - (ni + e*)) for all c G C \ {z}, (1) 

or cquivalcntly, 112 = z (ni + e,-) G CL(C) with z G C, thus ri2-ff = (ni + 
ei)H. Furthermore, the special case ofc = 0GC\{z}of Equation Q] implies that 
w_£r(ri2) < w^f(ni + e^). Therefore all condition of Definition HTl are verified, i.e. 
z = ni + n 2 + e, G L(C). 

Conversely, if z G L(C), then z can be rewritten as z = ni + ri2 + e.; where 

(1) m, n 2 G CL(C), 

(2) i £ supp(ni), 

(3) w H (m + ej) > w ff (n 2 ), 

(4) ffn 2 =%+e,). 

Now (1) and (2) gives that ni + e^ G X(D(0)), whereas (1), (3) and (4) clearly 
force that CL(ni + e.;) = ri2, i.e. wn(n-2) < wj/(ni + e, + c) for all c G C, or 
equivalently, nj + e, G D(z). Therefore X(D{0)) n D(z) ^ 0. □ 

Corollary 30. Let C be a binary code then Z m i n (C) C L(C), for any minimal test- 
set Z m - ln (C) obtained from Z(C). 

Proof. Let Z m i n (C) be a minimal test-set of C obtained from Z(C), then every 
z G Z milL (C) satisfies that ^(0(0)) n D(z) ^ 0. Thus, by Theorem [2^1 we obtained 
the required result. □ 

Remark 31. Algorithm [T8l gives the set of leader codewords L(C) of a binary code 
C. Furthermore any minimal test-set Z min is a subset of L(C), thus after performing 
redundancy elimination to L(C), also a minimal test-set iJ m i n can be obtained. 

Remark 32. There is a built in function CLBC2 in GBLA_LC which performs 
Algorithm [18] to compute L(C) and a variant of the same function which computes 
L\C). 

Example 33. We use the same code of Example [JS] Algorithm [T5] returns L(C) 
and L^C), in this case we got that both set coincides. We describe below the set of 
leader codewords with 14 elements of the given binary code C. 
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e 3 
e 2 
ei 



e 4 
e 3 
e 3 



e 7 
e 6 



68, 
67, 

e 7 , 



e 2 
ei 
ei 



e 4 
e 4 
e 2 



e 6 
e 5 



68, 

e6, 



e 4 + e 6 + e 7 + e 9 + ei , e 3 + e 6 + e 8 + e 9 + ei , 
L(C) = L X (C) = < e 2 + e 7 H 

ei + e 5 + e 6 + e 7 + e 8 + e 9 + eio, 
ei + e 3 + e 4 + e 5 + e 6 + e 9 + ei , 
6i + e 2 + e 4 + e 5 + e 7 + e 9 + eio, 
ei + e 2 + e 3 + e 5 + e 8 + e 9 + e i0 



Note that the only nonzero codeword of C that is missing in L(C) is the codeword 
y = (1,1,1,1,1,1,1,1,0,0) of weight 8. This result is consistent with the fact that 
the covering radius of C is p(C) = 3, as shown in Example 1151 and the statement of 
Theorem [24] where we proved that the weight of a leader codeword is always less or 
equal to 2p(C) + 1 = 7. 



(0,0,0,0,0,0,0,0,0,0), (1,0,0,0,1,1,1,1,1,1), (0,1,0,0,0,0,1,1,1,1), 

(1,1,0,0,1,1,0,0,0,0), (0,0,1,0,0,1,0,1,1,1), (1,0,1,0,1,0,1,0,0,0), 

(0,1,1,0,0,1,1,0,0,0), (1,1,1,0,1,0,0,1,1,1), (0,0,0,1,0,1,1,0,1,1), 

(1,0,0,1,1,0,0,1,0,0), (0,1,0,1,0,1,0,1,0,0), (1,1,0,1,1,0,1,0,1,1), 

(0,0,1,1,0,0,1,1,0,0), (1,0,1,1,1,1,0,0,1,1), (0,1,1,1,0,0,0,0,1,1), 



c = < 



(1, 1, 1, 1, 1, 1, 1, 1, 0, 0) 



In the following table we presents the computation results of a binary Golay 
code and a binary BCH code. 





[23, 12] Golay code 


[21, 12] BCH code 


Codewords (2 k ) 


4096 


4096 


Cosets (2"- fe ) 


2048 


512 


Leader codewords (|L(C)|) 


253 


478 


1^)1 


253 


549 



Table 1. Computation of the number of codewords, the number of cosets, the number of leader 
codewords and the cardinality of |L (C)| of the [23,12,7] binary Golay code and the [21,12,5] 
binary BCH code. 



Therefore we show an example where the subsets L(C) and L 1 (C) agree (this is 
not a surprise since the Golay code is a perfect code) and an example where the 
set L : (C) is smaller than L(C). Also note that both codes have the same number 
of codewords but the Golay code has four times the number of cosets of the BCH 
code. By contrast, the number of leader codewords is less in the Golay code. 
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5. Complexity Analysis 

For a detailed complexity analysis and some useful considerations from the compu- 
tational point of view we refer the reader to [8]. The main difference with respect 
to the previous works is that we not only provides a set of canonical forms AT but 
also the set of all coset leaders. The next theorem states an upper bound for the 
number of iterations that the Algorithm CLBC will perform. 

Theorem 34. The CLBC algorithm computes the set of coset leaders of a given 
binary code C of length n after at most n|CL(C)| iterations. 

Proof. Notice that by looking how the Algorithm [T] is constructed, the number of 
iterations is exactly the size of List. Moreover note that we can write List as the 
following set 

List = {w + e.j | w e CL(C) and i G {1, . . . , n}}. 
Therefore it is clear that the size of List is bounded by n|CL(C)|. □ 

Remark 35. 

(1) We can proceed analogously to the previous proof for estimating the required 
memory space which is O (n|CL(C)|). In the best case O (|CL(C)|) of memory 
space is needed so the CLBC algorithm is near the optimal case when consid- 
ering memory requirements. 

(2) The CLBC algorithm generates at most n|CL(C)| words from F£ to compute 
the set of all coset leaders. Therefore the proposed algorithm has near-optimal 
performance and significantly reduced complexity. 

(3) By construction the CLBC2 algorithm has the same time complexity as the 
CLBC algorithm. The advantage of computing the set of leader codewords is 
that it aid in solving similar problems as the function Matphi does but with a 
structure which is considerately smaller. 

6. Conclusion 

In this paper we have formulated the Algorithm CLBC which turns out to be 
quite efficient for computing all coset leaders of a binary code from the memory 
requirements and the computational complexity point of view although, as it is 
expected, the complexity of the algorithm is exponential (in the number of check 
positions). 

The difference of CLBC with its predecessors (see [3] and the references therein) 
relies also in the computation of all coset leaders instead of a set of representative 
leaders for the cosets, however, it supports the computation of the function Matphi. 
Note that the computation of Matphi is not necessary at all for the main goal of 
computing the coset leaders, we just have kept this resource in the algorithm because 
this structure provides some computational advantages (see j6!8l7j ) and, although 
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the algorithm will be clearly faster without computing Matphi, the nature of the 
computational complexity and space complexity will remain in the same order of 
magnitude. 

The Algorithm CLBC can be transformed into Algorithm CLBC2 in order to 
compute the set of leader codewords, which turn out to be a set of zeroncighbors 
that is a test set. The set of leader codewords allows to compute all coset leaders 
corresponding to the coset of a received word. We have shown also that the set of 
leader codewords contains any minimal test set. 

Unfortunately, a generalization to non-binary linear codes is not straigthforward 
from this work. The main reason seems to be that the solution based on Hamming 
weight compatible ordcrings will not continue being possible since the error vector 
ordering that is used in the general approach [8] is not a total degree compatible 
ordering, although it allows to set up the computational environment in order to 
compute Grobner representations for linear codes. 
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